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(57) Abstract 



i i i 
l 



Method of determining the stroke volume of the human heart from the pulse-type blood-stream pressure signal derived 
from the aorta and consisting in each case of a systolic and a diastolic period. Said method comprises calculating the flow q(t) 
from the pressure p(t) and integrating the flow over the systolic period, the aorta being regarded as a transmission line supple- 
mented with a windkessel compliance and the pressure/volume relationship in the aorta as an arctangent relationship. The pres- 
sure represented by the pressure signal - for each stroke prior to the calculation of the flow - is linearised by means of the arctan- 
gent relationship, for fixed aorta length, of the aorta pressure/cross section, the values, associated with the final diastolic 
pressure, of the characteristic impedance of the transmission line and windkessel compliance being adhered to or is used without 
modification. The characteristic impedance of the transmission line and the windkessel compliance are continuously adapted to 
the pressure of the windkessel compliance for each stroke concerned. Additionally, to calculate the flow, the flow through the pe- 
ripheral resistance is added. 
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Method of determining the stroke volume and the cardiac output of the 
human heart 

5 

The invention relates to a method of determining the stroke volume 
and the cardiac output of the human heart from the pulse-type blood- 
stream pressure signal derived from the aorta and consisting in each case 
of a systolic and diastolic period. It is known in practice to calculate 
* 10 the stroke volume and the cardiac output with the aid of a thermodilution 
determination by injecting cold liquid into the blood stream and 
measuring it downstream at regular intervals. 

In another known method, a pulse contour method or prescription is 
used to determine the stroke volume V s and the cardiac output or heart 
15 minute volume Q from a pressure signal p(t) measured in the human aorta. 
In this connection, the stroke volume is the volume of blood ejected by 
the heart in one contraction or stroke. A typical value is, for example, 
70 cm 3 , but this quantity can vary from stroke to stroke. The cardiac 
output is the volume of blood which the heart pumps in a unit time of one 
20 minute. In this connection, a typical value is 5 1/min, and this volume 
can typically vary between 2 and 30 1/min. 

In the past relatively simple relationships have been assumed 
between the measured aorta pressure signal p(t) and the mean aorta flow 
q(t). In this connection, the starting point is the pressure signal since 
25 said pressure can be measured relatively easily and well, but the flow 
cannot. In practice it has been found that such a relationship or "model" 
is much more complicated than was first assumed. This has emerged, in 
particular, in hospitals, where this so-called pulse contour calculation 
has resulted in frustration in the case of seriously ill patients. 
30 In practice, usually a so-called "Windkessel" or air-receiver model 

is assumed, the aorta being conceived as a single compliance (figure lb). 
In this case a windkessel - a container partially filled with liquid and 
with gas - which can absorb the surges in flow and buffer them in the gas 
bubble and which is incorporated downstream of a pulsatile pump is meant. 
35 The volume of blood ejected by the heart pump is largely received in the 
aorta or windkessel and partially flows away through the peripheral 
vascular regions of the various organs which branch off from the aorta. 
In the period in which the heart pump does not eject a volume of blood, 
i.e. the diastolic period T d , the outward flow from the windkessel 
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continues to be fed to the peripheral vascular regions. The aorta does 
not then contain any gas but the aorta wall is elastic and this elastic 
vascular wall fulfils the same function as the windkessel. 

This concept contains various factors which are neglected. Firstly, 
5 the aorta is much longer than wide. In fact, the pressure wave generated 
by the heart requires a time of 0.1 to 0.3 sec. to reach the end of the 
aorta and then approximately the same time to return to the heart, 
whereas the entire expulsion period of the heart (the systolic period or 
systole, T s ) lasts only 0.2 to 0.4 sec. At the beginning of the output 
• 10 flow, only a small portion of the compliance is therefore available for 
buffering the outward flow. However, still more important is the fact 
that, after two times the propagation or transit time of the pressure 
wave, twice the capacity is available, albeit with a doubling of the 
initial pressure wave amplitude. In addition, a second neglect which is 
15 generally made is that the windkessel is linear, that is to say a 

doubling of the stored stroke volume is accompanied by a doubling of the 
pressure. In reality it has been known for a number of years from the 
work of G.J. Langewouters et al. in J. Biomechanics 17. pages 425-435 
(1984) that the relationship between pressure p and volume V is strongly 
20 nonlinear in accordance with an arctangent function (figure 3a). In 

addition, said arctangent function is dependent on the age and sex of the 
patient. Finally, the drainage from the aorta to the peripheral vascular 
zones is not concentrated at one location but distributed along the 
aorta. That is to say, it takes a little time before the increased 
25 pressure at the start of the aorta has reached the various branching 
points to the peripheral vascular zones and the drainage actually 
increases due to the increased aorta pressure. 

To summarise, the neglects therefore amount to the fact that the 
aorta does not behave linearly and that there is a time factor due to a 
30 travelling wave along the aorta. 

In practice, however, another model is also assumed, the so-called 
transmission line model (figure la). In this approximation, the aorta is 
conceived as a homogeneous elastic tube of "adequate" length filled with 
liquid (blood) , with drainage to the peripheral vascular zones 
35 concentrated at the end. Such a tube has two characteristic properties. A 
pressure wave generated by the heart at the start of the tube travels at 
finite velocity v p to the end of the tube, is partially reflected there 
and travels back to the start of the tube. By the time the reflected wave 
arrives there, the expulsion phase, the systole, of the heart is already 
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over, the aorta valve is already closed and the heart no longer has any 
trouble from the reflected wave since the latter is, after all, held back 
by the closed valve. The "adequate length" is that length of tube which, 
given the propagation velocity of the wave, ensures that the reflected 
5 pressure wave does not return too soon. The propagation velocity v p is 
primarily determined by the area A of the cross section of the aorta, the 
compliance C of the aorta per unit of length and the density p of the 
blood in the tube as follows: v p = /A/pC . 

In addition to the propagation velocity, the second characterising 
• 10 property is the characteristic impedance Z 0 . This impedance reflects the 
ratio between the amplitude of the pressure wave and that of the 
accompanying flow wave in the aorta and is also affected by the cross 
sectional area A and the compliance C. In fact, the characteristic 
impedance is given by the formula: 
15 Z 0 = /p/(AC f ). 

The fact that the ratio between pressure p(t) and flow q(t) is fixed and 
given by the formula: q(t) = p(t)/Z Q means that the pressure wave and the 
flow wave at the start of the aorta are identical in shape at least till 
the instant that the reflected pressure (and flow) wave has (have) 
20 returned. 

This model also contains various approximations which, although 
they have less serious consequences for the stroke volume calculation, 
make the latter still too inaccurate under certain circumstances. 
Firstly, it is also true of this model that the aorta does not behave 

25 linearly in relation to the relationship between pressure p and volume V 
and therefore, with constant aorta length, as regards the cross section A 
= V/l. Since the compliance C per unit of length is the derivative of 
the area of the cross section with respect to pressure, at the prevailing 
pressure p c : C = [dA/dp]p Q# A will indeed increase with increasing 

30 pressure, but at the same time C will decrease because the aorta becomes 
increasingly less distensible with increasing pressure in accordance with 
the arctangent function mentioned (figure 3b). In a particular pressure 
region, these two factors will compensate for each other in the formula 
for Z Q and the latter will be virtually constant. At increasingly higher 

35 pressures, the decrease in C takes place more quickly than the increase 
in A, and the characteristic impedance Z c will start to slowly increase 
in accordance with the root of the inverse product of A and C . The 
nonlinearity therefore plays a role but this is less serious than in the 
windkessel model. A second approximation is that the aorta is assumed to 
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be homogeneous. In reality, the aorta is neither homogeneous in cross 
section nor in distensibility of the vascular wall. The aorta cross 
section decreases towards the periphery. This is initially compensated 
for because the arterial branches of the aorta start to contribute to the 
5 total cross section as soon as the pressure wave and flow wave reach the 
branches. In addition, the vascular wall is. however, increasingly less 
distensible towards the periphery. This manifests itself in a decrease in 
compliance C 1 . Since both A and C 1 in the formula for the characteristic 
impedance decrease over the aorta towards the periphery, Z 0 will 

10 increase. This means that there is no question of a single reflection at 
the peripheral resistance at the end of the aorta tube, as assumed in the 
model, but that there are distributed reflections along the entire aorta, 
initially small but considerably increasing later in the systole. With a 
short systolic period, the heart only experiences the (low) impedance at 

15 the start of the aorta while contracting, but with a longer-lasting 

systole it also experiences the (higher) peripheral impedance. This will 
mean that the instantaneously expelled quantity of blood will become 
increasingly smaller as the systole advances since it is gradually curbed 
by the increasing reflected impedance. 

20 A third approximation is the following. After the blood has been 

forced into the aorta due to an increase in the pressure, the pressure in 
the subsequent diastolic period will, of course, not be zero immediately, 
but the aorta will slowly empty through the peripheral resistance. As a 
consequence of this, the subsequent stroke will be forced out against a 

25 somewhat increased pressure and the two pressures will in fact be 

superimposed on each other. The diastolic pressure will always increase 
further with every stroke until equilibrium has been reached between 
systolic inward flow and diastolic outward flow. This is comparable to a 
windkessel function. 

30 Finally, it is also true that the drainage is not concentrated at 

the end of the aorta but is distributed over the length thereof. 

To summarise, the neglects are a consequence of the fact that the 
aorta has nonlinear properties and cannot be conceived as homogeneous. 
Both the above model concepts are not adequate if an attempt is 

35 made to calculate the aorta flow pattern from the aorta pressure as a 
function of time during a systole using the models. The curve found with 
the transmission line model only resembles an actual flow curve to any 
extent under the conditions of a short-duration systole. The windkessel 
model does not actually even permit a calculation of a flow curve. In 
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practice, for both models only the stroke volume of the heart, that is to 
say the flow curve integrated over one systole, is therefore calculated 
from the pressure curve. 

The object of the invention is to eliminate the abovementioned 
5 problems and to provide a method for the said determination of the stroke 
volume and cardiac output which takes account both of the windkessel 
properties and the transmission-line properties. This is achieved in a 
method, of the type mentioned at the outset, by calculating the flow q{t) 
from the pressure p(t) and integrating the flow over the systolic period, 
• 10 the aorta being regarded as a transmission line supplemented with a 

windkessel compliance and the pressure/volume relationship in the aorta 
as an arctangent relationship. In this case, the nonlinear aorta 
properties are important and are involved in this method, while it is 
advantageous for the verification of the correctness of the steps of the 
15 method if the flow curve can be continuously calculated and displayed 
during a systole. The integral of the flow curve then yields the stroke 
volume . 

The invention will be explained in greater detail on the basis of 
an embodiment, reference being made to the drawings, wherein: 
20 figures la and lb show the two known pulse contour models, viz. the 

windkessel model and the uniform transmission line model; 

figures 2a and 2b show the pulse contour model according to the 
invention based on the transmission line model supplemented by the 
windkessel function, respectively without and with peripheral drainage; 
25 figure 3a shows a graph with stylised pressure/cross section and 

pressure/compliance per unit length relationship for a 40 year old male 
aorta; 

figure 3b shows a diagrammatic representation of the calculation of 
a linearised pressure using the arterial final diastolic pressure as 
30 initial value; 

figures 4a and 4b show, respectively, an input pressure signal p(t) 
and a calculated flow curve q(t) calculated with the aid of a pulse 
contour model; and 

figure 5 shows a number of graphs of the pulse contour (line) and 
35 thermodilution (o) cardiac output of six open-heart operation patients 
plotted against the number of the thermodilution determination. 

The invention assumes the transmission-line model, in particular, 
the transmission-line property that pressure wave and flow wave are 
related by a fixed factor, the characteristic impedance Z c , at least 
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before reflections return to the heart from the periphery. This model 
lacks the gradual curbing of the outward flow from the heart by early 
reflections and the diastolic windkessel action. According to the 
invention, these two shortcomings are eliminated by means of adding one 
5 component, an aorta windkessel compliance C wf to the transmission-line 
model (see figure 2a). After all, in the systole, the flow is determined 
in the first instance by Z 0 and is then gradually curbed during the 
ejection phase by the pressure increase of p w (t) which is built up in the 
windkessel receiver by the inflow q(t). 

10 A further neglect was that the peripheral drainage was conceived as 

concentrated at the end of the transmission line. As a result of now 
associating, according to the invention, the equivalent peripheral 
resistance R p with the windkessel compliance C w , the latter is connected 
to an average (albeit low-pass -filtered) pressure and consequently 

15 represents reality better. See figure 2b. For in this model, the outward 

flow q R (t) is subtracted from the inflow q(t) so that the curbing of the 
p 

inflow takes place somewhat less rapidly. 

The improved model according to the invention contains linear time- 
invariant components. In reality, we are concerned with a strong 

20 nonlinear aorta whose pressure/cross section relationship varies in 

accordance with an arctangent. According to the invention two versions 
are now possible. In the first version, the pressure curve recorded in 
the method is first linearised and then applied to a linear pulse contour 
model. In the second version, the pressure measured in the method is 

25 regarded as given, but the pulse contour model is provided with pressure- 
dependent components. 

The linearised pressure curve model: 

In this first version it is assumed that the pressure curve arises 
because a certain volume of blood is pumped by the heart into the aorta, 

30 and this is followed by a rise in pressure. The magnitude of the rise in 
pressure depends on the pressure/volume relationship of the aorta. At a 
high pressure, a low volume increase will result in a high pressure 
increase, whereas at a low pressure, the same volume increase results in 
a smaller pressure increase (arctangent relationship). However, it is 

35 possible to calculate, by means of the pressure/volume relationship 

assumed to be known, the pressure curve which would have existed if the 
aorta had had a linear pressure/volume relationship. The procedure can be 
indicated graphically in a simple manner (see figure 3b). For every 
stroke, the start of the expulsion phase in the pressure curve and the 
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associated final diastolic pressure p d are determined. At said pressure, 
the nonlinear aorta pressure/ volume relationship or, assuming a fixed 
aorta length, pressure/cross section relationship is linearised by means 
of the tangent to the curve and the equivalent linear pressure is 
5 calculated for every pressure sample taken on the pressure curve. 

Assuming an initial pressure p d of 70 mmHg, in this example of a pressure 
sample of 100 mmHg, the equivalent linear pressure becomes equal to 
approximately 92 mmHg. This pressure is then applied to the pulse contour 
model which, for characteristic impedance and windkessel compliance 
* 10 calculated from the applicable pressure/cross section relationship, has 
the values associated with the final diastolic pressure for the stroke 
concerned . 

The nonlinear model: 

In this second version, it is assumed that the pressure curve 
15 arises as a result of the outward flow of the blood into a nonlinear 
aorta, given all its properties and not only its windkessel properties, 
as above. A linear correction is not then applied to the pressure curve 
but is accepted as it is. However, during the expulsion phase of the 
stroke concerned* the value of the characteristic impedance and of the 
20 windkessel compliance varies according to the pressure value such as it 
exists at any instant and is calculated from the nonlinear arctangent 
pressure/cross section relationship according to the above formulae. At 
any instant when a pressure sample is taken, the associated 
characteristic impedance and windkessel compliance, and the effect of 
25 these altered values on the flow are calculated. 

The use of the arctangent relationship: 

According to Langewouters et al., the relationship between the area 
of the aorta cross section and the aorta pressure can be specified using 
the formula: 

30 A(p) = AJ0.5 + (1/K.)tan^{(p-P 0 )/Pi}3 (D 

where A(p) is the cross section as a function of the pressure p, A ro is 
the limit cross section at very high pressure, and p Q and Pj are, 
respectively, the pressure at the point of inflection of the relationship 
and the slope of the relationship at the point of inflection. For a 
35 typical curve see figure 3a- 

The compliance C* per unit of length as a function of the pressure 
follows from this as a result of differentiation with respect to p as: 

C'(p) = [dA/dp] p = ^/{l+fp-pjVP! 2 } (2) 
where C m = A m /ff p = the maximum compliance at p = p 0 . The values of the 



WO 92/12669 PCT/NL92/00022 

8 

parameters A„, p c and p x are known from the work of Langewouters et al. In 
particular, these parameters depend on the sex and age of the person. The 
value for A m has, in addition, a dispersion of approximately 20%. It is 
not therefore possible to determine the values of the components in the 
5 pulse contour model in an absolute sense without more exact calibration 
unless a standard deviation of approximately 20% is regarded as 
acceptable . 

Z 0 and v p can readily be calculated with the aid of A and C found 
in this way. To calculate the total windkessel compliance of the aorta, 
• 10 still one more parameter is required, namely the effective length of the 
aorta. The windkessel compliance is then the compliance per unit of 
length multiplied by the effective length: C w = C 1 . -A- good value for for 
adults is 60 cm. 

The value of the peripheral resistance. 
15 The pulse contour model of figure 2b includes a component 

representing the total peripheral resistance R p of all the vascular 
drainage beds in parallel. The value of the latter is initially unknown. 
After the stroke volume V s has been calculated for the first time and the 
heartbeat rate /, the number of heartbeats per minute, has been found 
20 from the sum of the systolic and diastolic time periods T = T s + T d , if 
/ * 60/T, the cardiac output can be calculated as the product of the 
heartbeat rate and the stroke volume for the stroke concerned: 

Q = /V 6 . (3) 
The average aorta pressure p ro is also known as the integral of the 
2p pressure over a time period divided by the time period: 

P. « l/T)./p(t)dt CD 
T 

The value of the peripheral resistance then readily follows as: 

30 Rp - P»/0 (5) 

Since the peripheral resistance only varies slowly, its value, calculated 
for the present stroke, can be used for the flow calculation of the 
subsequent stroke. For the first stroke, a reasonable initial estimate 

35 will have to be taken. Since the value of the flow in the systole is not 
strongly dependent on the value of the peripheral resistance, a rapid 
convergence is obtained in this way. 

Calculation of the flow using the. pulse contour model. 

The flow q at any instant can be calculated using a simulation. In 

i|0 this connection, equations are set up which describe the behaviour of the 
model of figure 2b. The aorta flow follows from the equation: 
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q(t) = 1/Z e ).{(p(t) - P w (t)}. 

The flow through the peripheral resistance follows from: q R (t) = P w /R p ; 
and the volume in the air-receiver compliance C w follows from: 
v v(t) = V 0 + / {(q(t) - q R (t)}dt, where V G is the initial volume in C w . 
5 Finally, the pressure in the air-receiver follows from the 
volume/pressure relationship p w (t) = /(V w (t)). 

As stated earlier, the stroke volume V s can be found from the flow curve 
by integrating the flow during the systole, the heartbeat rate / as the 
quotient of 60 and the time period of every stroke, the cardiac output Q 

10 as the product of stroke volume V s and heartbeat rate and the 

peripheral resistance R p as the low-pass -filtered quotient of the average 
pressure p m and the cardiac output Q as indicated in formulae 3, 4 and 5. 

A measured pressure curve and the flow curve calculated from it are 
shown in figure 4 for a 62 year old male patient. Before every stroke, 

15 the beginning of the rising edge of the pressure curve in figure 4a is 
determined using a pattern recognition algorithm. At this instant, the 
diastolic pressure p d on the curve is recorded, and with the latter, the 
correct values for the characteristic impedance and the windkessel 
compliance are found using a function generator and introduced into the 

20 model. The pressure is linearised using the arctangent pressure/ cross 
section relationship and is applied to the model. The correct values for 
the parameters which determine the arctangent pressure/cross section 
relationship follow from regression equations according to Langewouters 
et al. and are as follows for women and men, respectively, having an age 
25 n: 

for women: for men: 

K ■ 4.12 A m = 5.62 

Po = 72 - 0.89n P c = 76 - 0.89n 

p, = 57 - OMn Pi = 57 - OMn 

30 where A,, is in cm 2 , the pressures p. is in mmHg and the age n is in years. 

Obviously, these calculations can advantageously be carried 
out by a digital processor. For this purpose, it is necessary to replace 
the continuous pressure curve by a series of samples with regular spacing 
intervals and to repeat the calculation for every sample taken. The flow 
35 curve therefore becomes available at these time instants. A preferred 
sample interval is 10 ms. A shorter interval of 5 ffi s is possible if a 
fast computer is used but only gives a marginal improvement in the 
precision. The integration of the flow curve in order to calculate the 
stroke volume is replaced by a summation of the flow samples in the 
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systolic period. 

The precision of the new method. 

The new method has been evaluated in six patients who were 
undergoing an open-heart operation. The reference value for the cardiac 
5 output was measured at some time instants during the operation by 

repeating a thermodilution (TO) determination at least four times after 
injections of cold liquid distributed with precise regularity over the 
respiratory cycle. A single thermodilution determination would be too 
inaccurate. The particular feature is that the four determinations were 
- 10 injected so as to be spread regularly over the respiratory cycle, as a 
result of which an averaging of the four values usually means an 
„ improvement in the precision of the estimate of the actual cardiac output 
by a factor of four. This improvement is only achieved with mechanical 
insufflation and not, therefore, if the patient breathes spontaneously. 
15 Spontaneous breathing always takes place during the first determination, 
i.e. at "s", as marked in figure 5. and during all the subsequent 
determinations, the patient is anaesthetised and he is artificially 
insufflated. 

The cardiac output calculated by the novel pulse contour method was 

20 recorded at the same time for comparison. The results are shown 
separately for each of the patients in figure 5- The number of 
simultaneous determinations, the mean error, the error standard deviation 
(SD) and the error span are shown for each patient number 8910 to 8912 
and 8915 to 8917 in table 1 and averaged for the group. Note that the 

25 sixth patient has a larger error for three determinations. These were 

carried out after a period of total heart-long bypass while, although the 
heart was already beating and pumping blood , the pump of the heart-lung 
machine was still assisting at approximately one litre per minute. This 
flow was in fact measured by the thermodilution method but not by the 

30 pulse contour method, resulting in this discrepancy. 

For comparison, the cardiac output was also calculated by the 
earlier method according to Wesseling et al. from US 3.841,313. which 
method was known as the most accurate in the literature. The mean error 
standard deviation was in this case 12*, against 8% for the novel method. 

35 The novel method therefore has a variance of less than half that of the 
old. 

Roughly speaking, all the mean errors - shown in the last two 
rows - with the new method are more than a factor of /2 lower than with 
the older method. 
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Table 1 



pat nr 


TD 


PC 




Error old 




Error new 






old 


new 


mean 


SD 


span 


mean 


SD 


span 


8910 


4.67 


5.21 


4.53 


0.54 


0.43 


1.02 


0.14 


0.38 


1.12 


K>y±X 




4 8Q 


t • yy 


0.54 


0.38 


1 02 


0.44 


0.37 


0.99 


8912 


4.59 


5-31 


5.41 


0.73 


0.45 


1.24 


0.08 


0.38 


1.37 


8915 


4.10 


4.49 


3.80 


0.38 


0.65 


1.81 


0.30 


0.26 


0.60 


8916 


4.69 


4.98 


5.13 


0.29 


0.45 


1.44 


0.43 


0.33 


0.89 


8917 


4.29 


4.33 


4.03 


0.04 


0.88 


2.97 


0.26 


0.50 


1.50 


mean 


4.63 


4.87 


4.50 


0.24 


0.54 


I.58 


0.13 


0.37 


1.08 


SD 


0.46 


0.39 


0.52 


0.45 


0.19 


0.74 


0.30 


0.08 


0.33 



During the intensive monitoring of a patient in the operating theatre 
or intensive-care ward, it is only rarely possible to record the arterial 
pressure in the aorta. For reasons of convenience and safety for the 

20 patient, the pressure in the a. radialis or brachialis or a finger 
pressure will normally be recorded non-invasively. This pressure is 
distorted compared with the aorta pressure, that is to say the shape of 
the pressure curve has changed. This is due to the fact that the pressure 
curve generated in the aorta propagates through the arterial system to 

25 the more peripherally situated upper arm, pulse or finger and that this 
transmission path has certain properties, the most important of which are 
its narrow-band nature and its low-frequency resonance. It is possible to 
measure these properties in a large number of patients. In this 
connection, the inventor has found that the transmission properties have 

30 a relatively small spread around a mean value. In that case it is 

possible to correct the resonance properties by an an ti -resonance in a 
filter, as a result of which it is indeed possible to accurately 
determine the pulse contour stroke volume of the left-hand half of the 
heart. To a good approximation, a simple second-order filter can be used 

35 for this purpose. In practice, a digital, IIR (infinite impulse response) 
second- or third-order filter is used which corrects the mean resonance 
precisely up to a frequency of 8 Hz. 

In this case, on the basis of a series of input signal samples x n , a 
series of output signal samples y n is calculated using the formula: 
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x„. een reeks van uitgangssignaal monsters y n berekend met de formule: 
H 

y„ = c 0 k„ * E d k y n . k 
k=l 

5 where c and d are constant filter coefficients. For the second-order 
filter, the following values are used: 
c 0 = 0.077 
d, = l.W* 
d 2 = -0.561 
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Claims 

1. Method of determining the stroke volume of the human heart 
from the pulse-type blood-stream pressure signal derived from the aorta 
and consisting in each case of a systolic and a diastolic period, 
5 comprising calculating the flow q(t) from the pressure p(t) and 

integrating the flow over the systolic period, the aorta being regarded 
as a transmission line supplemented with a windkessel compliance and the 
pressure/volume relationship in the aorta as an arctangent relationship. 

10 2. Method according to claim l f wherein the pressure, represented 

by the pressure signal, is linearised prior to the calculation of the 
flow using the arctangent relationship, with fixed aorta length, of the 
aorta pressure/cross section in accordance with 
A(p) = A m [0.5 * (VnJ.tan-^fp-pJ/p,}], where 
15 A(p) = the cross section as a function of the pressure p, 
Ag, = the limit cross section at very high pressure, 
p D = the pressure of the point of inflection of the arctangent 
relationship, and 

Pj = is the slope of the arctangent relationship at the point of 
20 inflection. 



3. Method according to claim 2, wherein the pressure linearised 
for each stroke is used to calculate the flow, and wherein, for each 
stroke concerned, the values, associated with the final diastolic 
25 pressure, of characteristic impedance of the transmission line and 
windkessel compliance are adhered to. 



4. Method according to claim 1, wherein the pressure, represented 
by the pressure signal, for each stroke is used to calculate the flow, 
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and wherein, for each stroke concerned, the characteristic impedance of 
the transmission line and the windkessel compliance are continuously 
adapted to the pressure of the windkessel compliance. 



5 5. Method according to claim 3. wherein the flow q(t) is 

calculated according to . 

q(t) = (l/Z 0 ).{p(t)-p w (t)}, where Z G = characteristic impedance 
= /p/AC), and p w » the compliance pressure, and 
p = density of blood, 
10 C'(p) « C ro /{l+(p-po) 2 /Pi 2 }. where 

C„ = ^/npj = the maximum compliance at p « p D . 

6. Method according to claim 1, wherein the cardiac output Q is 
calculated as the product of the stroke volume and the heart rate. 

15 

7. Method according to claim 6, wherein, additionally, to 
calculate the flow, the flow through the peripheral resistance R f is 
added which follows from 

q R (t) = p v /Rp. where R p = p m /Q, where p m = mean pressure = 
p 

20 (1/T) Jp(t)dt. 



8. Method according to one of the preceding claims, wherein, if 
the blood pressure is measured in the arteria radialis or brachialis or 
non invasively in the finger, a filter operation is performed on the 
25 pressure signal in order to correct the resonance distortion in the 
pressure curve caused by the transmission. 
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